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ABSTRACT 

We show that asymmetries in total intensity and Hnear polarization between the radio jets 
and counter-jets in two lobed Fanaroff-Riley Class I (FRI) radio galaxies, B2 0206+35 
(UGC 165 1) and B2 0755+37 (NGC 2484), can be accounted for if these jets are intrinsically 
symmetrical, with decelerating relativistic outflows surrounded by mildly relativistic back- 
flows. Our interpretation is motivated by sensitive, well-resolved Very Large Array imaging 
which shows that both jets in both sources have a two-component structure transverse to their 
axes. Close to the jet axis, a centrally-darkened counter-jet lies opposite a centrally-brightened 
jet, but both are surrounded by broader collimated emission that is brighter on the counter- 
jet side. We have adapted our previous models of FRI jets as relativistic outflows to include 
an added component of symmetric backflow. We find that the observed radio emission, after 
subtracting contributions from the extended lobes, is well described by models in which de- 
celerating outflows with parameters similar to those derived for jets in plumed FRI sources 
are surrounded by backflows containing predominantly toroidal magnetic fields. These return 
to within a few kpc of the galaxies with velocities ^ 0.25c and radiate with a synchrotron 
spectral index a ^ 0.55. We discuss whether such backflow is to be expected in lobed FRI 
sources and suggest ways in which our hypothesis can be tested by further observations. 
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1 INTRODUCTION 



Relativistic jet outflows from radio galaxies are a primary mecha- 
nism for energy extraction from supermassive black holes in ac- 
tive galactic nuclei (AGN) and an important source of energy 
input to the intergalactic medium (IGM) in groups and clus- 
ters (e.g. iMcNamara & Nulsenll2007l and references therein). We 
are studying relativistic jet kinematics and dynamics in nearby 
low -luminosi ty radio g alaxies with Fanaroff-Riley Class I (FRI - 
iFanaroff &Rilevlll974l) morphology for which we have obtained 
radio imaging and polarimetry at high angular resolution trans- 
verse to the jets as well as along their lengths. We have developed 
procedures for deriving three-dimensional variations of intrinsic 
jet parameters - velocity field, emissivity and magnetic-field or- 
dering - from an an alysis of systemati c asymmetries between the 
jets and counter-jets (iLaing & Bridle 2002a [lCanvin & Laind2004l : 
(Canvin et aD[2003;|Laing. et al.i i2006a). We compare the observed 
asymmetries in images of total intensity, degree of linear polariza- 
tion and apparent magnetic field direction with the predicted effects 
of relativistic aberration on synchrotron emission from particles in 
partially-ordered magnetic fields in model outflows and deduce the 
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distributions of intrinsic properties within the jets. We have found 
that a generic property of the jet outflows in FRI radio galaxies 
is that they decelerate from relativistic speeds {(3 — v/c ^ 0.^ - 
0.9) near the AGN to subrelativistic speeds a few kiloparsecs away, 
and that the outflows are systematically faster on-axis than at their 
edges. 



It is critical for such an analysis to distinguish patterns of 
asymmetry in the jets produced by relativistic aberration from any 
that are intrinsic to the outflows or which result from interactions 
between the outflows and anisotropic environments, e.g. from pres- 
sure gradients or winds in the IGM. One asymmetry in FRI radio 
jets that has proven instructive in some sources and problematic 
in others is the systematic difference between transverse intensity 
profiles in the brighter jets and wealcer counter-jets when observed 
at high sensitivity and angular resolution. 



This difference correlates with indicators of the orientation of 
the jets to t he line of sight. A statistical study of FR I jets in the B2 
sample by Laing et al. (1999) found that the ratio of jet to counter- 
jet FWHM measured by Gaussian fitting at the same distance from 
the nucleus on both sides is strongly anticorrelated with the aver- 
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age jet/counter-jet brightness ratio and with the ratio of cor^B to 
extended flux density. 

This anticorrelation is quaUtatively as expected for intrinsi- 
cally symmetrical relativistic outflows which are faster on-axis 
than at their edges. In this case, relativistic aberration makes the 
transverse brightness profiles of the approaching, hence appar- 
ently brighter, jet more centrally peaked than those of the reced- 
ing counter-jet. Gaussian fitting to the jet and counter-jet FWHM 
then yields smaller values of the width for the apparently brighter 
jets even if the (slower moving) outer boundaries of the jets appear 
identical on both sides of the AGN. 

The ampli t ude o f the effect found in the B2 source sample 
by iLaing et aP (Il999h is, however, surprisingly large. Modelling 
of the anticorrelati on requires that th e velocity /3on-axis ~ 0.7 
and /Sedge ~ 0.1 jLaing et al.l ll999h in order to reproduce the 
spread of width ratios. Two lines of argument suggest that such 
large velocity ratios are not typical of the FR I population. Firstly, 
the ratio /3edge //^on-axis required to explain the effect is quanti- 
tatively inconsistent with the brightness and polarization distribu- 
tions in four of the five individual FR I sources we have modelled 
(Eaing & Bridlel2002al : ICanvin & Laing 2004l: ICanvin et al.l2005l - 

the exception is 3C 296: lLaing et al 2006a) . Secondly, the smallest 

values of /3edge //^on-axis arc required only to generate the unusu- 
ally small values of jet/counter-jet width ratio ^ 0.6 in a few mem- 
bers of the B2 source sample with particularly high jet/counter-jet 
brightness ratios, whose jets are tho ught to be highly inclined to the 
plane of the sky (iLaing et al.ll 19991) . 

Thus far, our results would be consistent with the idea that 
all FRI jets are symmetrical outflows, but that only a few have 
very large transverse velocity gradients. Even this hypothesis fails 
for two of the B2 sample members, B2 0206+35 and B2 0755-^37 
jLaing et al.ll20TlF[ These sources are unusual in that the lower 
isophotes of their brighter jets also appear narrower than those 
of the counter-jets at the same distance from the AGN in images 
of moderate resolution and sensitivity (e.g. iBondi et al ]|200P) - 
even though the jets clearly exhibit the basal asymmetries associ- 
ated with symmetrical decelerating relativistic outflows. Apparent 
width asymmetry in the fainter jet emission cannot generally be ex- 
plained by relativistic effects alone if the jets are both symmetrical 
and purely outflowing On the other hand, if the asymmetry is at- 
tributed to intrin sic or environment al differences on the two sides 
of the AGN (e.g. 'Bondi et al.|[200Q ^ there should be no systematic 
trend for the wider j et to be on the re ceding side as it is in the (albeit 
small) sample of iLaing et al.l (Il999l) . 

In this paper, we explore an alternative explanation for the 
transverse brightness profile asymmetries of the jets and counter- 
jets in 0206-1-35 and 0755-1-37. This work was motivated by new 
deep imaging of these sources showing: (a) that their counter-jets 
have minima in their emission profiles with the same widths as the 
main jets at similar distances from the nucleus and (b) that the main 
jets are surrounded by faint emission resembling the broader outer 
emission in the counter-jets (ILaing et al .1 12011, and Section [2l2l 
below). The new imaging data lead us to model the jets in these 
sources as intrinsically symmetrical outflows near the jet axis sur- 
rounded by broader features from bacMowing material. If back- 



flow in the broader features can be approximately symmetrical and 
mildly relativistic, then aberration can make its emission appear 
slightly brighter on the counter-jet side, producing differences in 
isophotal width between the jets similar to those observed. 

Backflow is a reasonable hypothesis a priori for FRI sources 
like 0206-1-35 and 0755-1-37 whose jets appear to propagate within 
well-defined lobes. It has been an acknowledged ingredient of mod- 
els of lobed FR I I sources since the fi rst attempts to simulate their 
hydrodynamics (lNormanetal.lll982h . FRI sources cannot form 
lobes without similar deflection of jet material and iLaing et aP 
(120111) showed that FR I lobes resemble those of FR II sources in 
many respects. If FR I jets are much lighter t han their surroundings 
and initially fast (e.g. iLaing & Bridlel2002bh . we should not be sur- 
prised if some large-scale post-jet flow in FRI lobes is marginally 
relativistic. We also note that mildly relativistic backflow extends 
almost all the way back to the centre of the host galaxy in sim- 
ulations of relativistic FRI jets with initial dynamical flow pa- 
rameters matching those deduced from our observations of 3C 31 
and realistic pressure and density profiles for the surrounding IGM 
(Laing & Bridle 2002b; Perucho & Marti 2007). 

In this paper, we show that a fully symmetrical model in which 
a decelerating axisymmetric outflow is surrounded by a slower (but 
still slightly relativistic) backflow is quantitatively consistent with 
the detailed brightness and polarization distributions of the jets and 
counter-jets in 0206-1-35 and 0755-1-37. It is not obvious a priori that 
conditions needed to produce symmetrical backflow are likely to 
be realised in lobed FRI radio galaxies. Nevertheless, our results 
suggest that mildly relativistic backflow contributes significantly to 
the observed jet vs counter-jet width relationships and we suggest 
ways in which this (perhaps unexpected) ingredient of FR I source 
structure could be investigated further. 

In Section [21 we summarize the optical and large-scale ra- 
dio properties of the sources and discuss the additional image pro- 
cessing required to separate jet and lobe emission. Section [3] de- 
scribes our modelling procedure and Section |4]gives a comparison 
between models and data. The model parameters are presented in 
Section |5] A brief discussion is given in Section [6] Section|7] sum- 
marizes our conclusions and suggests further work. Finally, Ap- 
pendix |A1 demonstrates that a toroidally-magnetized outflow can, 
in special circumstances, produce jet/counter-jet sidedness ratios 
significantly less than unity. 

We adopt a concordance cosmology with Hubble constant, Hq 
= 70kms"^ Mpc"\ Qa 0.7 and Om 0.3. 



2 OBSERVATIONS AND IMAGES 

2.1 The sources: optical data and large-scale radio structures 

The galaxy identifications, redshifts and linear scales for the two 
sources studied here are given i n Table [T] T heir radio structures 
have been described in detail by iLaing et aP (12011.) . from which 
the images in Fig. [T] are taken. 



^ The 'core' is defined as an unresolved component coincident with the 
AGN. The core/extended flux-density ratio is a statistical indicator of ori- 
entation. 

^ From now on we drop the B2. 

^ We discuss a special magnetic-field configuration for which this is not 
the case in Appendix IaI 



2.2 Images 



Table |2] summarizes the relevant parameters of the high- 
resolution sub-images which we model or use for spectral anal- 
ysis in this paper (details of the observations and data reduction 
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Table 1. Names, redshifts, linear scales and associated references for the 
sources in this paper. 



0.5 1.0 1.5 2.0 2.5 



Name 


Galaxy 


Redshift 


Scale 


Reference 




name 




kpc 










arcsec"-*^ 




0206+35 


UGC1651 


0.03773 


0.748 


1 


0755+37 


NGC2484 


0.04284 


0.845 


2 



References: (l) lMilleret 011120021) : (2) lFalco et alii 19991) . 
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Figure 1. Grey-scale images of the sources (iLaing etal.l l201lh . The 
boxes mark the areas shown in later plots and the grey-scale ranges, in 
mJy beam" ^ , are indicated by the labelled wedges, (a) 0206+35 at 4.9 GHz, 
1.2 arcsec FWHM. (b) 0755+37 at 4.9 GHz, 1.3 arcsec FWHM. 



are given by i Laing et al .1120 111) . The E- vector position angles of 
linear polarization at 4.860 GHz have been corrected for Faraday 
rotation using mu ltifrequency imaging dGuidetti et aT]l201lLl2012l : 
iLaing etalJl2oTll) and residual depolarization is predicted to be 
negligible at this frequency. The areas plotted in later figures are 
outlined on Fig.[T] 

Fig. [2 shows rotated sub-images. On the right-hand side of 
each panel, we have drawn a single contour to outline the brightest 




0.5 1.0 1.5 2.0 2.5 




Figure 2. False-colour images of total intensity for 0206+35 and 0755+37 
over the areas outlined in Fig. [T] On the right-hand side of each panel, we 
have plotted a single contour to outline the brightest emission of the main 
(brighter) jet. On the left-hand side, this contour (rotated through 180°) 
is plotted on the counter-jet emission, (a) 0206+35 at 0.35-arcsec FWHM 
resolution, (b) 0755+37 at 1.3-arcsec FWHM resolution. 



part of the main jet. On the left-hand side, this contour is rotated 
through 180° and plotted on the counter-jet emission. This diagram 
emphasizes the points made earlier that the minima in the counter- 
jet emission have roughly the same widths as the main jets and that 
the main jets are in turn surrounded by fainter emission. 

In order to model jets that appear superimposed on lobes, we 
must try to separate the two emission components in all Stokes 
parameters. There is no unique way to do this when their spectra 
and intensities vary independently across the field of view. Any ap- 
proach to isolating jet emission in a lobed FRI source therefore 
entails some simplifying assumption about the variations in inten- 
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Table 2. Parameters of the sub-images used for modelling and spectral anal- 
ysis. Col. 1: source name; col. 2: observing frequency (an asterisk indicates 
that the image was used for modelling); col 3: resolution (FWHM); col. 4: 
rms off-source noise level in /; col.5: average noise level in Q and U\ col. 6: 
sub-image position angle; col. 7: sub-image sizes parallel and perpendicular 
to the jet axis. 



Source u Res rms Rot Size 

GHz arc- /iJyb"-*^ deg arcsec^ 







sec 


(71 


crp 






0206-^35 


1.425 


1.20 


19 




-41.0 


22 X 20 


0206-^35 


4.860 


1.20 


12 




-41.0 


22 X 20 


0206-^35 


4.860 


0.35* 


7.2 


7.1 


-41.0 


22 X 20 


0755-^37 


1.425 


1.30 


20 




158.5 


66 X 66 


0755+37 


4.860 


1.30* 


7.8 


7.9 


158.5 


66 X 66 


0755-h37 


4.860 


0.40* 


8.0 


7.1 


158.5 


20 X 16 



sity / or in spectral inde^fl a of the lobes or jets over the region 
to be modelled. We have attempted to separate the jets and lobes 
for these sources in a way that optimizes the resolution and signal- 
to-noise of the jet emission while letting us check for systematic 
errors resulting from the assumptions made while doing the lobe- 
jet separation, as follows. 

One approach to separating jet and lobe emission observed at 
two frequencies is based on their systematic spectral differences: 
the jets have characteristic spectral indices close to a = 0.55, 
whereas the lobes have a > 0.8 near the centres of the sources 
(iLaing et aP l201lh . If the spectral index of the lobe emission 
close to the jet is reasonably cons tant, we can use a variant of 
the 'spect ral tomofflap hy' metho d (Katz-Stone & Rudnickl ll997l : 
iKatz-Ston e et al. 1999; Laing et al. 2006b) by assuming that what 
is observed can be described as the sum of two components: a jet 
and a lobe with constant spectral indices and a\, respectively. 
The brightnesses observed at a given point at two frequencies 
and vi are then: 

We can scale and subtract the two brightness distributions to esti- 
mate the jet brightness at the modelling frequency uq: 

Once we know a\, the method can also be applied to Stokes Q 
and U provided that we correct the images at both frequencies for 
Faraday rotation before subtraction, and that depolarization is neg- 
ligible (as is the case for these sources). Note that the spectral index 
of the jets, aj, must be both constant and known in order to scale 
the result correctly. 

In practice, we estimated the lobe spectral index for vq = 
4.860 GHz and z^i = 1.425 GHz by performing the subtraction for 
various trial values of ai and selecting that which minimized the 
residual lobe emission in jet-free regions. 

Spectral subtraction can remove even rather complicated lobe 
emission if the spectral index is constant, but it has two serious 
flaws for our purposes: (a) the signal-to-noise ratio of the corrected 

^ We define spectral index a in the sense I{h') oc 



Table 3. Interpolation parameters for lobe subtraction. Col. 1 : source name; 
col 2: resolution (FWHM); col. 3: background region distances from jet 
axis; col. 4: width of boxcar smoothing function parallel to the axis. 



Source 


FWHM 


Background 


Smooth 




(arcsec) 


(arcsec) 


(arcsec) 


0206-^35 


0.35 


9-10 


1.0 


0755-^37 


1.30 


30-45 


3.0 



image is lower than that of the deep high-frequency image alone 
and (b) our highest-resolution data for 0206-1-35 and 0755-1-37 are 
only at one frequency. 

The alternative of spatial subtraction assumes that the lobe 
intensity varies only slowly across the jet. This approach can be 
best applied at high angular resolution where the lobe brightness 
is low and the spatial variation of jet emission is clearest. To sep- 
arate the two types of emission spatially in /, Q and [/, we define 
two background regions parallel to the jet axis and just outside the 
maximum transverse extent of the jet as estimated from spectral- 
index images, i.e. using both the intensity and spectral properties 
of the jet emission to guide our choice of the background regions. 
We then smooth the background brightness distributions parallel to 
the jet axis with a boxcar function to improve their signal-to-noise 
ratio and interpolate linearly between them under the jetjf] We refer 
to this approach as generating 'interpolated images'. 

For 0206-1-35 and 0755-1-37 we first used spectral subtraction 
to verify the total extent of the jet emission and to set appropriate 
reference regions for interpolation, then constructed interpolated 
images for the final modelling. 

In Figs[3]and|4l we show the results of both subtraction meth- 
ods for the two sources. Figs|3ja) and Ha) show the images at the 
resolution used for modelling before subtraction. We found best- 
fitting lobe spectral indices between 1.425 and 4.860 GHz of 0.90 
and 0.81 for 0206-^35 and 0755-^37, respectively. In Fig.|3lb), we 
show the spectral subtraction for 0206-1-35 at lower resolution. Al- 
though not useful for modelling, this image outlines the total extent 
of the flatter- spectrum emission associated with the jets. The spec- 
tral subtraction for 0755-1-37 at the lower of the two resolutions used 
for modelling, shown in Fig. |Hb), has little trace of residual lobe 
emission but low signal-to-noise. 

Guided by the spectral subtraction, we set the interpolation pa- 
rameters as in Table |3] and computed interpolated images at 1.425 
and 4.860 GHz, from which we in turn derived the spectral-index 
images shown in Figs^c) andUc). These are blanked on the er- 
ror in spectral index, as noted in the captions. We then estimated 
integrated spectral indices for the jets by summing the interpolated 
/ images over all pixels which are unblanked on the spectral-index 
images, excluding the cores. We found (aj) = 0.55 for 0206-1-35 
and 0.53 for 0755-1-37. We used these values to scale the spectral 
subtractions. Variations across the modelled regions are small, with 
0.50 < ttj < 0.62 in both sources. Finally, we show the 4.860- 
GHz interpolated images at the resolutions used for modelling in 
FigsHd) andHd). 

We are confident that the interpolated images represent the 
jet emission accurately in both sources. The lobe emission in 
0206-1-35 is quite faint at 0.35-arcsec FWHM resolution, and af- 
ter subtraction, the area around the jets appears devoid of resid- 

^ Higher-order interpolation works poorly for these brightness distribu- 
tions. 
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(d) Interpolated 







Figure 3. False-colour images showing the results of lobe subtraction for 
0206+35. The / intensity colour range (0-3 mJy beam~^) is the same for 
panels (a) and (d). (a) No subtraction at 0.35-arcsec resolution, (b) Subtrac- 
tion at 1.2-arcsec resolution assuming a constant spectral index for the lobe, 
(c) Spectral index distribution over the jet and counter-jet at 1.2-arcsec reso- 
lution after interpolated subtraction, blanked where da. > 0.03 (the colour 
range for spectral index is shown by the labelled wedge), (d) Subtraction by 
linear interpolation between background strips parallel to the jet axis. The 
resolution is 0.35 arcsec. 

ual emission in all Stokes parameters (e.g. Fig. |3}i). The lower- 
resolution (1.3 arcsec FWHM) image of 0755+37 proved to be 
more of a challenge, because the lobe emission is bright and ir- 
regular (Fig. |4t). The spectral subtraction gave a clean image of 
the jet with negligible background emission, but amplified noise 
(Fig. IIId). In contrast, interpolation (Fig. |4ji) failed to remove the 
small-scale lobe emission accurately but retained the full signal- 
to-noise ratio of our high-frequency images. Comparison of the 
two corrected / images showed that they are accurately consistent 
wherever / > 100 /xJy beam~^ . We therefore used the interpolated 
images for modelling (in which the faint residual lobe emission 
has low weight). Modelling the spectrally- subtracted image (and its 
counterparts in Q and U) gave consistent but less well constrained 
results. In the intensity and polarization profiles plotted below, we 
compare the results from both subtraction methods. 

At 0.4-arcsec FWHM resolution, used for modelling the inner 
jets of 0755+37, the lobe brightness is negligible and we did not 
attempt to subtract it. 



3 MODEL FITS 
3.1 Assumptions 

To model the jet emission, we make the following assumptions. 

(i) The jets are intrinsically symmetrical, axisymmetric and an- 
tiparallel. They can be treated, on average, as laminar, stationary 
flows. 

(ii) The radio emission is from relativistic particles with a 
power-law energy spectrum n{E) — noE~^'^^~^'^^ (a is the spec- 
tral index). We use the integrated values for the modelled regions 
after lobe subtraction: (a) = 0.55 for 0206+35 and 0.53 for 
0755+37. The corresponding maximum degree of polarization is 
po = (3a + 3)/ (3a + 5) = 0.70 in both cases and the variations 
of spectral index across the modelled regions are small enough to 
be ignored (Section I2l2l ). 



0.0 0.5 1.0 1.5 2.0 


0.4 0.5 0.6 0.7 0.8 
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(c) Spectral index 


A 




^^^^^^^ 


■:::3 


(b) Spectral 
subtraction 




(d) Interpolated 




1 




1 arcsec 



Figure 4. False-colour images showing results of lobe subtraction in 
0755+37 at 1.3-arcsec resolution. The / intensity colour range (0 - 
2.5 mJy beam~^) is the same for panels (a), (b) and (d). (a) No subtraction, 
(b) Subtraction assuming a constant spectral index for the lobe, as described 
in the text, (c) Spectral index distribution over the jet and counter-jet after 
interpolated lobe subtraction, blanked where a a. > 0.1 (the colour range 
for spectral index is indicated by the labelled wedge), (d) Subtraction by 
linear interpolation between background strips parallel to the jet axis. 

(iii) The magnetic field is tangled on small scales, but 
anisotropic. 

(iv) The effects of Faraday rotation on the observed emission 
are corrected completely. T his is an extreniely good approximation 
for 02 06+35 and 0755+37 (iGuidetti et al.ll201lLl2012l : lLaing et alJ 

[mil). 

3.2 Outline of method 

For a symmetrical, outflowing jet with velocity v = pc, emitting 
isotropically in the rest frame and inclined by an angle to the line 
of sight, a measurement of the observed jet/counter-jet intensity 
ratio 

= [(1 + /3 cos ^)/(l - /3 cos ^)]'+" 

does not allow us to determine the velocity and inclination sepa- 
rately. The key to our method is the use of linear polarization to 
break this degeneracy. The relation between the angles to the line 
of sight in the rest frame of the outflow, 0' and in the observed 
frame, 0, is: 

sin^j = [r(l -/3cos6')]"^sin6' (mainjet) 
sin 6*^] = [r(l + /3 cos e)]~^ sin (counter-jet) 

The emission in all three Stokes parameters depends on 0' , since the 
magnetic field is in general anisotropic. If the flow is significantly 
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relativistic, we effectively observe the two jets at different values of 
0' and can use the differences in polarization for the approaching 
and receding jets as an additional constraint to separate [3 and 0. 
For backflow, the argument is identical with the roles of jet and 

counter-jet interchanged. 

The principal ste ps in our method (Laing & Bridle 2002a; 
[Canvin & Laing. 20041 : ICanvin et al...2005 : Laing et al...2006a) are 
as follows. 

(i) Build a parameterized model of the geometry, the velocity 
field and the variations of emissivity (oc noB^^") and magnetic- 
field anisotropy in the rest frame of the emitting plasma. 

(ii) Calculate the observed-frame emission in /, Q and U, taking 
account of relativistic aberration and anisotropic emission in the 
rest frame. 

(iii) Integrate along the line of sight, normalize to the measured 
total flux density and convolve with the observing beam. 

(iv) Calculate and sum over the /, Q and U images. This is 
our measure of goodness of fit. 

(v) Optimize th e parameters using the downhill simplex method 



(V) Optimize th e parameters using tJ 
of Nelder & Mead jPress etal. Ill992h . 



We explored a wide range of starting simplexes in order to be sure 
of locating the global minimum in x^. 




3.3 Fitting functions 

The parameterized model that we fit to the VLA o bservations is a 
simpli fi ed version of those in our previous wo rk ("Laing & Bridle" 
2002a| : ICanvin & Lain3 l2004l : ICanvin et al.l |2005; Laing et al. 
2006ah . with the addition of a few extra terms to describe the back- 
flow. The functional forms are given explicitly in Table |4] A criti- 
cal discussion of fitting functions will be given elsewhere (Laing & 
Bridle, in preparation). 



3.3.1 Geometry 

We use coordinates (z, x) in a plane containing the jet axis, with z 
measured along the axis and x perpendicular to it. The jet is divided 
into SL flaring region, where the flow first expands and then recolli- 
mates, and a conical outer region, as sketched in Fig. \5\ The edge 
of the outflow is fully defined by the distance of the transition be- 
tween the two regions measured along the axis, ro, the radius, xo, 
and the opening angle of expansion in the outer region, ^o- Indi- 
vidual streamlines in the outer region are straight, so we can define 
a streamline index s = ^/^o, where ^ is the angle between the 
streamline and the axis, s ranges from on-axis to 1 at the edge 
of the outflow. The two coefficients a2(s) and 03(5) of the cubic 
expression for the streamline radius in the flaring region (Table |4]) 
are defined by the conditions that the streamline radius x{z) and 
its first derivative x\z) are continuous at the flaring-outer region 
boundary. We also define a distance coordinate r which is contin- 
uous along a given streamline from at the nucleus to ro at the 
flaring-outer region boundary and which thereafter increases as the 
distance from the boundary surface. The functional forms for r in 
the two regions are given in terms of z and s in Table lU 

For simplicity, the backflow is assumed to follow the same 
streamline family as the jet, extended away from the axis. The edge 
of the backflow in the outer region is defined by the radius Xh at 
the region boundary and the opening angle ^b- These are not inde- 
pendent: Xh/xo = sin^b/ sin^o- The backflow streamline index 
t ranges from at the backflow/outflow interface to 1 at the edge 



Figure 5. Sketch of the assumed geometry. The blue and green curves show 
the outer boundaries of the outflowing jet and backflow emission, respec- 
tively, Representative streamlines in the two parts of the flow are shown in 
red. The fiducial distances and angles are defined in Section |3".3. II 



of the outflow. The backflow streamline radii have the same func- 
tional form as their outflow equivalents with the coefficients a2 (t) 
and as (t) again defined by continuity at the region boundary. 

The assumed backflow geometry is ad hoc, but gives a reason- 
able match to the observed extent of the emission. 



3.3.2 Velocity 

The on-axis velocity profile in the outflow is divided into three 
parts: (a) constant with a high velocity close to the nucleus; (b) 
a linear decrease and (c) constant with a low velocity at large dis- 
tances. The velocity along any off-axis streamline is calculated us- 
ing the same expressions but with truncated Gaussian transverse 
profiles. The velocity profiles, given explicitly in Table |4] depend 
on two transition distances, rvi and r^o, the on-axis velocities /3i 
and /3o and the fractional edge velocities vi and vo (which are re- 
quired to be < 1). 

We experimented with several functional forms for the back- 
flow velocity. The most satisfactory has no dependence on r, but 
varies linearly with streamline index from /3b, in at the interface with 
the outflow to /3b, out at the outer edge of the backflow. 



3.3.3 Emissivity 

We write the proper emissivity as e/, where e is the emissivity in 
Stokes / for a magnetic field perpendicular to the line of sight and 
/ depends on the field geometry (defined in Section [3". 3. 41 below), 
e, to which we refer loosely as 'the emissivity', is a function only of 
the total rms magnetic-field strength and the normalizing constant 
of the radiating electron energy distribution. 

The on-axis emissivity profile in the outflow is also divided 
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Table 4. Coordinate definitions and functional forms for geometry, velocity, proper emissivity and magnetic-field ordering. 



Description 



Quantity Functional form 



Distance range 



Distance coordinate r 
(outflow and backflow) 

Outflow streamline index s 



Outflow radius 



Outflow velocity 



x{z^ s) 



(ro+A)cos^-A 
cos 4 

A = xo/sin^o - r-o = Xb/sin^b - '^o 



by continuity 

(2; - ro + xq/ sin^o) tan(^o5) 

/3i exp(s'^ In vi) 

ft exp(.2l„„i) (-1^) +^oexp(.2ln^o) {^^) 

l3o exp(s2 In vq) 



Outflow proper emissivity ^{r,s) gir exp(ln eis"^) 

I exp l^ln — 



(^eO-^) + eo(^-^el) 



Outflow (^2/^2) 1/2 ji^r) 



Outflow (Bf/^?)^/^ k{r) 



gor ^ou* exp(lneo52) 



Ji 

jl{rBO-r)+jo{r-rBl) 
rBO—rBl 

jo 

ki 

ki(rBO-r)-\-ko{r-rBi) 
rBO-rBl 

ko 



r < To 
r > ro 



r < ro 
r > ro 

r < ro 
r > ro 

r < rv-L 

Tvi <r < rvo 

r > 

r < rei 

rei <r < reo 

r > reO 

r < TBI 

TBI <r< rso 

r > rso 

r < rBi 

rBi <r< rso 

r > tbo 



Backflow streamline index t 



Backflow radius 



x{z^ t) 



Backflow velocity /3 {t) 

Backflow proper emissivity e(r, t) 



Backflow {^2/^2^1/2 ■ 
Backflow (Bf/Bf) 1/2 



by continuity 

«-^o)/«b-^o) 

a2{t)z'^ + a3{t)z^ 

- ro + xo/ sin^o) tan[^o + Kb - ^o)*l 

/5b, in + *(/5b,out - /5b, in) 



n\){r/ro) exp(lnebt2) 

Jh 
kh 



r < ro 
r > ro 

r < ro 
r > ro 



r < ri 
r >ri) 



into three regions, each with a power-law profile. The profile is 
allowed to be discontinuous at each of the region boundaries. Off- 
axis, the profile is multiplied by a truncated Gaussian function of 
the streamline index, with values at the jet edge which are constants 
in the inner and outer emissivity regions and vary linearly between 
them. The free parameters for the emissivity profiles are transition 
distances, Veo and rei, power-law indices Ein, E^id and ^out, gi 
and go, which measure the discontinuities at the region boundaries 
and edge emissivities ei and eo. Note that ei and eo may be > 1 
(in which case the jet is limb-brightened), = 1 (uniformly filled) or 
< 1 (centre-brightened). 

The backflow emissivity is assumed to be zero within a given 
distance and to have a power-law dependence on r with a single 
index and a truncated Gaussian dependence on streamline index t 
elsewhere. The fitted parameters are the index E\,, the fractional 
edge emissivity Cb, the inner distance rb and the emissivity ratio 



between outflow and backflow at the boundary between the flaring 
and outer regions, rih. 



3.3.4 Magnetic -fie Id structure 

We define the rms components of the magnetic field to be {Bfy^"^ 
(longitudinal, parallel to a streamline), (5^)^^^ (radial, orthogo- 
nal to the streamline and outwards from the jet axis) and {B^y^'^ 
(toroidal, orthogonal to the streamline in an azimuthal direction). 
The rms total field strength is 5 = {Bf + B^ + 5?)^^^ The 
magnetic-field structure is parameterized by the ratio of rms ra- 
diaytoroidal field, j = {B^^^^ / {B^y^'^ and the longitudi- 
nal/toroidal ratio k = {BfY^'^ /{Bt)^^'^ . For the outflow models in 
the present paper, these depend only on r, being constant close to 
and far from from the nucleus and varying linearly at intermediate 
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distances. The free parameters are the fiducial distances tbi and 
tbq and the field ratios at these distances, ji, jo, ki and /cq. 

For the backflow, we assume constant field ratios jb and kh. 

3.4 Modelling of individual sources 

We estimated the noise levels for each resolution and Stokes pa- 
rameter based on the deviations of the brightness distributions from 
those expected for axisymmetry, as follows. 

(i) Calculate Stokes parameters Q and [/ in a coordinate system 
with position angle along the jet axis. 

(ii) For / and Q, take the noise level to be 1/ times the rms 
difference between the image and a copy of itself reflected across 
the jet axis. 

(iii) For U, take the sum rather than the difference. 

These values can be substantially larger than the off-source rms, 
but include the effects of small-scale structure (which we do not 
attempt to model) and deconvolution errors. 

For 0206-1-35, we fit to images at the highest available res- 
olution, 0.35 arcsec FWHM, using different noise levels for the 
high-brightness emission close to the nucleus and the fainter re- 
gions farther out. For 0755-1-37, we fit to 0.4-arcsec FWHM images 
of the bright inner jets and 1.3-arcsec images elsewhere. Small re- 
gions around the cores were excluded from the fits, since we model 
only optically-thin emission. The model images given below in- 
clude point sources with the appropriate observed flux densities at 
the locations of the cores. 

The values of summed over all Stokes parameters and res- 
olutions were 8012 over 6696 independent points for 0206-1-35 and 
8022 over 5816 points for 0755-^37. 

The quoted uncertainties were also derived as in our earlier 
work by varying an individual parameter until increased by an 
amount corresponding to the formal 99 per cent confidence level, 
leaving the rest of the model unchanged. These values are crude 
(they neglect coupling between parameters), but in practice give a 
good impression of the range of reasonable models. As an addi- 
tional check, we also performed a series of optimizations at fixed 
values of and tabulate the range over which acceptable solutions 
could be found. 



4 MODEL-DATA COMPARISONS 
4.1 General 

In Figs[6]-[8]and[T0l-[T2l we show various comparisons between 
the observed and model images of the two sources. The images 
have been rotated by the angles given in Table |2] so that the main 
(approaching) jet points to the right and the core is either at the cen- 
tre or the left-hand edge of a plot. The types of plot are as follows. 

(i) False-colour images of total intensity. The angular scale is 
given on the accompanying profiles and the brightness range (in 
mJybeam"^) is indicated by the labelled wedges. 

(ii) Longitudinal profiles of total intensity. 

(iii) Images of jet/counter-jet sidedness ratio /j//cj derived by 
dividing the / image by a copy of itself rotated by 180°. These 
images are blanked (grey) where / < Scrj on either side of the core 
(Table [3). The contours show /j//cj = 1. Angular scales are again 
shown on the accompanying profiles. 

(iv) Longitudinal profiles of sidedness ratio. 



(v) Images of degree of polarization, p = P/^- These are 
blanked wherever / < 5cr/. The angular scale is given on the 
accompanying profiles and the range is indica ted by the labelled 
wedg es, p has been corrected for Ricean bias (iWardle & Kronbera 
|l97i. 

(vi) Profiles of p along the jet axis. 

(vii) Vectors with lengths proportional to p and directions along 
the apparent magnetic field, superposed on false-colour images of 
/. The angular and vector scales are indicated by labelled bars. 

(viii) Averaged transverse profiles of total intensity, /, sidedness 
ratio /j / /cj , and Q/I over selected regions where the brightness 
and polarization distributions vary slowly with distance from the 
nucleus. Stokes Q is defined in a coordinate system with its axis 
along the jet: Q/I > for an apparent magnetic field transverse 
to the axis; Q/I < for a longitudinal field. In the flaring region, 
these profiles were derived by averaging along radii from the nu- 
cleus, in which case they are plotted against angle from the jet axis. 
For the outer region, they are averages along lines parallel to the jet 
axis and are plotted against angular distance from the axis. In order 
to make a fair comparison, only pixels which were not blanked on 
the observed images were used in the averages. 

In general the fits are very good. We examine the correspondence 
between model and observed brightness distributions in detail in 
the next two sub- sections. 

4.2 0206+35 

We show images and longitudinal profiles of total intensity and sid- 
edness ratio in Fig. [6] and of degree and direction of polarization in 
Fig- El Averaged transverse profiles of /, /j//cj and Q// are given 
in FigM 

The model accurately reproduces the main features of the 
brightness and polarization distributions of 0206-1-35, including the 
following. 

(i) The main (approaching) jet has a bright base, with a peak at 
^2 arcsec from the nucleus (Figs|6t - c). 

(ii) The peak sidedness ratio of /j / /cj ~ 37 is at a distance of 
^^0.6 arcsec from the nucleus (Figs|6}i - f), close to the position 
of the flaring point as determin ed from high-resolution MERLIN 
observations jLaing et al.ll201lh . 

(iii) At low isophotes, the counter-jet appears wider than the 
main jet (Figs[6t and b). 

(iv) The counter-jet has a limb-brightened structure, which is 
brightest between 2.5 and 6 arcsec from the nucleus, whereas the 
main jet appears narrower and is centrally peaked (Figs |6^ and b; 
Figsit-f). 

(v) The longitudinal profile of degree of polarization shows the 
characteristic asymmetry we have noted in other FR I jets: the main 
jet has a polarization minimum at ^2.5 arcsec from the nucleus, 
corresponding to the transition between longitudinal and transverse 
apparent field, whereas the counter-jet shows a high degree of po- 
larization with a transverse apparent field, reaching an average of 
p 0.5 at 10 arcsec (Fig.|7t). 

(vi) There is a transition in the field direction between transverse 
on-axis and aligned with the jet boundaries at the edges on both 
sides of the nucleus. This is clear within 2 or 3 arcsec of the ridge 
line in the main and counter-jets, respectively (Fig.|7]and Figs [8] - 
o). The signal-to-noise ratio in the data is too low to determine the 
edge field direction accurately at larger distances, so discrepancies 
between observed and predicted Q/I transverse profiles should not 
be taken too seriously. 
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Figure 6. Comparison between the observed and modelled total intensities / and sidedness ratios /j//cj for 0206+35. (a) observed and (b) model false-colour 
images of /. (c) profiles of observed (full/red) and model (dashed/blue) / along the axis of the jet. (d) and (e) images of /j //cj . The white contours represent 
/j//cj = 1: outside the contours, /j//cj < 1- (f) profiles of observed (full/red) and model (dashed/blue) /j//cj along the jet axis. 



(vii) Close to the nucleus, the apparent field wraps around the 
edges of both jets, with a high degree of polarization, especially on 
the counter-jet side (Figs|7li and e). 



The main deficiency of the model is that it underpredicts the 
brightness of the counter-jet ^ 5 arcsec from the axis and over- 
predicts that of the main jet between 1.5 and 4 arcsec. These effects 
lead to a model sidedness ratio which is too high off-axis, although 
still significantly < 1 . This discrepancy is most obvious between 5 
and 7 arcsec from the nucleus (Figs[8j), e and h), but is restricted to 
regions where the brightness is < 200 /iJy beam~^. The model is 
also constrained to have monotonic deceleration in the outflow and 
velocity independent of distance from the nucleus in the backflow, 
so it cannot reproduce the increase in sidedness ratio between 8 and 



10 arcsec from the nucleus. The surface brightness is low at these 
distances so uncertainties in lobe subtraction may be significant. 

Fig. |9] shows the predicted brightness distributions for the 
outflowing and backflowing parts of the model separately. The 
former is similar to the pure outflow n i odels we have derived 
for other sources (iLaing &Bridlell2002at |c anvm &LaindE00i 
ICanvin et aLll2005l : iLaing et al.ll2006ain n the model of 0206+35, 
the limb-brightening of the counter-jet is due to a combination of 
outflow and backflow. In the outflow, the on-axis velocity remains 
high, so the edges of the outflowing counter-jet material appear rel- 
atively brighter because they suffer less Doppler dimming than the 
on-axis material. This effect is reinforced by emission from the 
backflow, which adds a thin shell of emission immediately sur- 
rounding the outflow. Most of the asymmetry is due to the outflow: 
the backflow is only slightly brighter on the counter-jet side. 
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Figure 7. Comparison between the observed and modelled linear polarization of 0206+35. (a) and (b) colour images of degree of polarization p = P/I in the 
range - 0.7, as indicated by the labelled wedge. Blanked areas are grey, (a) observed; (b) model, (c) profiles of observed (full/red) and model (dashed/blue) 
p along the axis of the jet. (d) and (e) vectors with lengths proportional to p and directions along the apparent magnetic field superimposed on colour images 
of /. (d) observed, (e) model. 



4.3 0755+37 

We compare model and observed total intensity images and profiles 
for 0755+37 in Fig.[T0l the corresponding polarization comparisons 
are shown in Fig.[TT]and Fig. [T2l gives averaged transverse profiles 
of /, /j / /cj and Q/I. Note that the fainter emission is affected by 
imperfect lobe subtraction, as discussed in Section 1231 This is par- 
ticularly obvious at large distances from the jet axis in images of 
ratios such as /j / /cj and p. The following features of the brightness 
and polarization distributions are reproduced. 

(i) The main jet has a brightness peak at 1 .3 arcsec from the core 
(Figs[ToJ - k). Farther out, the profile declines rapidly with distance. 

(ii) There is a rapidly-expanding, triangular region of roughly 
uniform brightness at the base of the counter-jet fFigslTOt and b). 

(iii) The jet base structure is initially very asymmetric, with a 
peak sidedness ratio ?^80 at 1.9 arcsec from the core, decreasing 



rapidly with distance to reach an asymptotic value at 15 arcsec 
(Figs[T0t-h). 

(iv) At faint brightness levels, the counter-jet appears signifi- 
cantly wider than the main jet, with a large opening angle (FigsfTOk 
and b). 

(v) The counter-jet brightness profiles are more flat- topped or 
edge-brightened than those of the main jet at most distances from 
the core (FigsfTok and b and Figs[T2k - h). 

(vi) A prominent arc of emission crosses the counter-jet at 
^26 arcsec from the nucleus (FigsfToh and b). 

(vii) There is also a bar of emission crossing the counter-jet at 
^12 arcsec from the nucleus (FigsfTOh, b and d). 

(viii) The profiles of degree and direction of polarization along 
the axis show the same characteristic asymmetry seen in 0206+35 
and other FRI jets. There is a change in apparent field direction 
at ^5 arcsec from the nucleus in the approaching jet, but not in 
the counter-jet, whose apparent magnetic field is always transverse 
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Figure 8. Transverse profiles of total intensity, /, jet/counter-jet sidedness ratio, Ij/hj and Q/I for 0206+35. The data have been averaged parallel to the 
jet axis over three ranges of distance from the nucleus: 3 - 5, 5 - 7 and 7 - 9 arcsec, as indicated in the captions (Section ITU . Full/red line: observations; 
dashed/blue line: model. Q/I > and Q/I < correspond to transverse and longitudinal apparent field, respectively. 
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(b) 0206+35 backflow 



ever, exaggerated by the look-up table in Figs [Toja) and (b) and is 
more accurately represented by the profile in Fig. [T2la). 

Fig. [13] shows the outflow and backflow components of the 
model intensity distribution. As for 0206+35, the outflow appears 
similar to that in other FR I radio galaxies, but the backflow is rela- 
tively stronger in 0755-1-37. The prominent curved arc crossing the 
counter-jet ^26 arcsec from the nucleus is modelled as the projec- 
tion of the inner edge of the backflow at r = rb. This is roughly 
elliptical in shape, with an axial ratio of sec — 1.22 and there is 
good correspondence between model and data. As mentioned ear- 
lier, the fit to the bar crossing the counter-jet closer to the nucleus 
is less successful. In the model, this is the other half of the pro- 
jected inner edge of the backflow, so there is no freedom to adjust 
its location or curvature to match the observed feature more closely. 
A similar problem afflicts the main jet: the projection of the inner 
edge of the backflow appears slightly too bright, causing the excess 
off-axis emission close to the nucleus. 



5 arcsec 



5 DERIVED PARAMETERS 

The best-fitting parameters for our models of 0206-1-35 
0755-^37 are Usted in Tables [5] (outflow) and [6] (backflow). 



and 



Figure 9. Predicted brightness distributions for the outflowing and back- 
flowing parts of the model for 0206+35. (a) outflow; (b) backflow. 

(Figs[TT^ and h, Figs[T2Ji - t). The degree of polarization in the 
counter-jet rises monotonically with distance from the nucleus, 
reaching large values {p ^ 0.5) far from the nucleus (Figs [TTt 
-c). 

(ix) The degree of polarization in the main jet base is low, and 
the apparent field is longitudinal (Figs[TTJi - f, i and j). 

(x) There are minima in the degree of polarization on either side 
of the axis in both jets, corresponding to the transition between 
transverse and longitudinal apparent field (FigsITT^, b, g, h:[T2ln - 
t). 

(xi) There is a region of high polarization with a circumferential 
magnetic field around the base of the counter-jet (Figs fTT^ and h). 

(xii) Determination of the observed polarization in the faint re- 
gions far from the axis is complicated by imperfect subtraction of 
lobe emission, but the apparent field is primarily parallel to the 
edges of both jets (Figs[T2jn - 1). 

Features which are not fit well by the model are as follows. 

(i) The observed brightness distribution of the bright main jet 
base is slightly more centre-brightened than the model and the ob- 
served degree of polarization is higher than predicted at its edges 
(Figs[T0l,j;[IIli, e, i,j). 

(ii) The observed transverse total-intensity profiles are signifi- 
cantly more limb-brightened than the model in some places, and in 
particular between 18 and 21 arcsec on both sides of the nucleus 
fFigs[T2b and g). 

(iii) The inner bar crossing the counter-jet is both straighter and 
slightly farther from the nucleus in the observed image (?^13 arcsec 
compared with ?^11 arcsec for the model; Figs [Tok, b, d). The fit 
may be affected by the limb-brightening in this region, however. 

(iv) As in 0206+35, the off-axis brightness of the main jet is 
slightly overestimated close to the nucleus. The difference is, how- 



5.1 Geometry 

Both sources are fairly close to the line of sight, as expected from 
their high jet/counter-jet sidedness ratios and bright cores. We de- 
rive = 40° for 0206+35 and 35° for 0755+37. The outflow ge- 
ometries are typical of those we have determined for other FRI 
jets, with the boundaries between flaring and outer regions at 5.3 
and 13.9 kpc from the nucleus for 0206+35 and 0755+37, respec- 
tively. The corresponding half-opening angles in the outer regions 
are 3? 9 and 7? 4. 

In 0206+35, the backflow has a half-opening angle of 11° in 
the outer region and its emission extends back into the flaring re- 
gion, with a cut-off at rb = 2.7 kpc. For 0755+37, on the other 
hand, the backflow emission is truncated within the outer region 
(rb = 23 kpc), where its half-opening angle is 16°. 



5.2 Velocity 

Velocity images derived from our model fits are shown in Fig. [141 
The initial velocities of both outflow components are similar 
(/3i = 0.86 for 0206+35 and 0.88 for 0755+37) and the asso- 
ciated transverse velocity profiles are close to uniform. 0206+35 
shows little on-axis deceleration, reaching an asymptotic velocity 
/3o = 0.68 after 4 kpc. Its transverse velocity profile evolves much 
more, and the fractional edge velocity is 0.04 at large dista nces. 
In bot h these respects, the source resembles 3C296 (L aing et al.l 
l2006ah . 0755+37, on the other hand, appears to decelerate rapidly, 
to /3o = 0.25 by 18.5 kpc, with a fractional edge velocity of 0.26. 
This estimate should be treated with caution since the emission 
in the outer counter-jet is dominated by the backflow component, 
making it difficult to assess the intensity or polarization of the out- 
flow there. 

The backflow velocities increase away from the source axis, 
from 15 = 0.05 to 0.20 for 0206+35 and from 0.25 to 0.35 for 
0755+37. 
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Figure 10. Comparison between the observed and modelled total intensities / and sidedness ratios /j//cj for 0755+37. (a) and (b) colour images of /. 
(a) observed; (b) model, (c) and (d) profiles of observed and model / along the axis of the jet. (e) and (f) images of /j //cj . The white contours represent 
/j//cj = 1: outside the contours, /j//cj < 1- (g) and (h) profiles of observed and model /j//cj along the jet axis. The resolution for panels (a) - (h) is 
1.3 arcsec FWHM. (i) and (j) colour images of / in the range 0-5 mJy beam~^ for the base of the main jet. (i) observed, (j) model, (k) profile of observed 
and model / along the jet axis. The resolution for panels (i) - (k) is 0.4 arcsec FWHM. In the profile plots, the full and dotted (red) lines are from observed 
images with lobe subtraction by interpolation (as for the colour plots) and spectrum, respectively. The dashed/blue line is the model. 
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Figure 11. Comparison between the observed and modelled linear polarization for 0755+37 at resolutions of 1.3 and 0.4 arcsec FWHM. (a) and (b) colour 
images ofp = P/I in the range - 0.7 at 1.3 arcsec FWHM. (a) observed; (b) model, (c) profiles of observed (full/red) and model (dashed/blue) p along the 
axis of the jet. Only the profile of p derived from interpolated images is plotted; the equivalent for spectral subtraction is very noisy, (d) - (f): as (a) - (c) but 
for the main jet only at 0.4 arcsec FWHM. (g) and (h): vectors with lengths proportional to p and directions along the apparent magnetic field superimposed 
on colour images of /. The resolution is 1.3 arcsec FWHM and the vector scale is indicated by the labelled bar. (g) observed, (h) model, (i) and (j): as (g) and 
(h), but for the main jet at 0.4 arcsec FWHM. 
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Figure 12. Transverse profiles of total intensity, /, jet/counter-jet sidedness ratio, /j//cj and Q// for 0755+37. The data have been averaged along radii from 
the nucleus from 8.1 - 12.6arcsec and from 12.6 - 17.1 arcsec and parallel to the jet axis from 18-21 arcsec and 24 - SOarcsec, as indicated in the captions 
(Section im . Full and dotted (red) lines both represent observations, with lobe subtraction by interpolation and spectral methods, respectively. Dashed/blue 
lines show the model. Q/I > and Q// < correspond to transverse and longitudinal apparent field, respectively. 
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tative argument, the best-fitting transverse emissivity profiles are 
higher at the edges than on-axis. This effect is shght in 0206+35, 
where the profile is consistent with a uniformly-filled cylinder ev- 
erywhere. In 0755-1-37, however, limb-brightening is required over 
much of the outer region (Fig.fTSb). As noted in Section B31 the ob- 
served transverse intensity profiles in this source are significantly 
more limb-brightened than the model predicts, suggesting that there 
is a narrow enhancement in emissivity at the boundary between 
the outflow and backflow. The functional form we assume for the 
transverse variation of emissivity does not allow for such narrow 
features. 

The backflow emissivity decreases with distance at similar 
rates in the two sources (oc r~^"^^ in 0206-1-35 and oc r~^"^^ 
in 0755-^37). It is centre-brightened in 0206-^35 (eb = 0.02) but 
closer to uniform in 0755-1-37 (cb = 0.79). 





5.4 Field Ordering 

The fractional components of magnetic field, {B^Y^'^/B 
(toroidal), {B^f'^jB (longitudinal), and {Blf'^jB (radial) are 
plotted in Fig.[T6l 

In both sources, the field close to the nucleus in the out- 
flow is close to isotropic, with the longitudinal component just ex- 
ceeding the other two. At larger distances, the toroidal component 
dominates, with significant longitudinal and radial contributions in 
0206-1-35 and 0755-1-37, respectively. As for velocity and emissivity, 
the field components in the outer parts of 0755-1-37 may have larger 
systematic errors because of the dominance of backflow emission. 

The field in the backflow is toroidally dominated in both 
sources, with non-negligible radial components in both cases and 
some longitudinal field in 0206-1-35. 



Figure 13. Predicted brightness distributions for the outflowing and back- 
flowing parts of the model for 0755+37. (a) outflow; (b) backflow. 



5.3 Emissivity 

Model images of n^^B^^^ (proportional to the emissivity function 
e) are shown in Fig. [15] 

The model outflow components again show properties very 
similar to those in other FRI jets. The locations of the flar- 
ing points (0.82 and 1.55 kpc from the nucleus for 0206+35 and 
0755+37, respectively) are well determined and consistent with 
higher-resolution observations (Laing et al. 2011). The emissivity 
variations in the faint and poorly resolved inner jets upstream of 
the flaring points are not well constrained. In the flaring and outer 
regions, the gradient of the emissivity profile flattens with distance 
in both sources, as is usual in FRI jets. 0755+37 requires a sudden 
decrease in emissivity with distance at r = reo whereas 0206+35 
does not. 

The observed limb-brightening in both sources shows side-to- 
side symmetry. This cannot result from a transverse velocity gra- 
dient in the sense we have inferred, which would lead to limb- 
brightening only in the counter-jet. In agreement with this quali- 



5.5 Backflow spectral index 

We can also constrain the spectral index of the radio emission from 
the backflows. The spectral indices at the edges of the jets, where 
the line of sight is mainly through backflow emission after the lobe 
subtraction, are much closer to those of the jets themselves than 
to the values elsewhere in the lobes. We can estimate the spectrum 
of the backflow emission directly from the images shown in Fig- 
ures [3lc) andlUc) or, more accurately, by integrating total intensity 
at 1.425 and 4.860 GHz over pixels which are unblanked in these 
images. The latter method gives mean spectral indices of 0.50 for 
0206+35 and 0.57 for 0755+37, compared with 0.55 and 0.53 for 
the sum of outflow and backflow emission. 



6 DISCUSSION 

6.1 Testing the hypothesis 

It is clear that the initial jet base asymmetries of most FRI Jets 
are produced by rela t ivistic aberration (Laing & B ridle 2 00 2al : 
ICanvin&La"in3i2004 ICanvin et all l2005l : iLainglTalJ l2006ah . If 
0206+35 and 0755+37 prove to be typical - in that counter-jets 
consistently appear wider than the main jets at a given isophote 
in lobed FRI sources whose jet base asymmetries are large - then 
the jet/counter-jet width asymmetry must also be correlated with 
jet orientation. The models presented in Section|4]show that mildly 
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Table 5. Model parameters which are common to outflow and backflow, or Table 6. Model parameters for backflow (Section |33] and Table |4). Col. 1: 

which apply only to the outflow (Section [33l and Table|4). Col. 1: parameter; parameter; col. 2: unit; cols 3 and 4: values for 0206+35 and 0755+37. 

col. 2: unit; cols 3 and 4: values for 0206+35 and 0755+37. The parameters 

are defined in Section [33] and fisted in Table|4] is the range of angles to 

the line of sight for which any acceptable solutions can be obtained. Variable 0206+35 0755+37 



Variable 



0206+35 



0755+37 



Geometry (common to outflow and backflow) 



Geometry 



deg 
kpc 



10.9 



+0.5 
0.5 
9 7+0.1 



15 
23.2 



+0.5 
-0.1 
+0.8 
0.7 






deg 


40.0+O-J 


34 8+°-'^ 






AO 


deg 


34-43 


32.5 - 37.5 






ro 


kpc 


5 o+O.l 




/5b,in 












/5b, out 








Outflow geometry 






Co 


deg 




7.4ig:? 






XQ 


kpc 


-I 09+0.02 


QQ+0.08 


rib 


xlOO 



rvo 

/31 

^0 
VI 

vo 



kpc 
kpc 



Velocity 



2+0.3 
^-0.3 
1+0.3 
-0.2 
0.86+0-08 



lA 
4.V 



-0.07 
3+0.09 
-0.05 
- + 0.05 
" -0.13 

^•^^-0.01 



0.68: 
0.95^ 



3.6 



+1.6 



18.5^ 



-1.5 

p; + 2.3 

0.88+Si 

-0.05 
1.00-0.06 

^•^"-0.11 



0.25^ 



^b 
eb 



Jh 



Velocity 



0.02 
0.20 



+0.03 
0.02 
+ 0.06 
0.07 



Emissivity 



9 0+0.2 
^•'^-0.2 

-^•""-0.07 



r) 9P'+0.04 

n op: + 0.05 
^•'^^-0.05 



0.094 



+ 0.000 
-0.010 

1.81+0-0^ 



0.79 



0.05 
+0.13 
0.14 



Field component ratios 



0.24^ 



^+0.08 
-0.07 

38+0-08 



0.38 
0.03 



+0.07 
0.07 
+0.15 
0.03 



Emissivity 



reO 
Eout 

ei 
eo 
91 
90 



TBI 

rBO 

31 

jo 

ki 

ko 



kpc 
kpc 



0.82 
2.04 



+0.02 
0.02 

+ 0.07 
0.06 

^ 3.1 

2 59+0-0^ 

^•"^^-0.08 



2.13 



+0.08 
-0.06 
1 9+0.6 

1 14+0.16 
-^•-^^-0.16 
1 7+0.8 
' -1.3 

1 n^+o-os 



1.55 
10.2 



+0.04 
-0.03 
+ 0.1 
0.3 

^ 2.4 

o 7f^+0.02 
'"-0.04 



1.16 



+0.05 
-0.09 

1.0+0-3 



2.2 
1.7 
0.52 



0.2 
+0.5 
0.3 
+0.5 
0.4 

+ 0.06 
0.03 



Field component ratios 



kpc 
kpc 



4.6 
1.50 
0.11 
1.36 
0.64 



< 1.4 

+0.5 
0.5 
+0.34 
0.22 

+ 0.13 
0.11 
+ 0.13 
0.13 
+0.05 
0.04 



Q Q + 2.8 

1 5 4+2-5 

^+0.13 
-0.09 



0.96^ 



0.44 



+0.12 

0.15 
+0.08 



1 1 t=;"r 
0-08io.08 



relativistic backflow offers a possible cause for such an orientation- 
dependent effect. 

There is an alternative explanation for /j//cj becoming <1 
in some parts of a source which also preserves the orientation- 
dependence of the effect. For the special case where the magnetic 
field is purely toroidal and the edge velocity is ^ cos ^, it is pos- 
sible for relativistic aberration to give an off-axis jet/counter-jet 
sidedness ratio < 1 even for a pure symmetrical outflow. We anal- 
yse this special case in Appendix [A] where we show that it is in- 
consistent with the polarization imaging of 0206+35 and 0755+37. 




Figure 14. The model values of velocity /3 in units of c in planes containing 
the jet axes. Positive and negative values of /3 denote outflow and backflow, 
respectively, (a) 0206+35, (b) 0755+37. 
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(b) 0755+37 




Figure 15. The model values of loginoB^^^) in planes containing the jet 
axes (no and B are in SI units), (a) 0206+35, (b) 0755+37. 



The mechanism inevitably produces degrees of polarization close 
to the theoretical maximum of po ~ 0.7 with a transverse appar- 
ent field. It is therefore unlikely to be important in the majority of 
observed jets but it may be relevant in a few objects like 3C296 
(Appendix lAt. 

If the jets are intrinsically symmetrical, then the backflow hy- 
pothesis remains the most plausible explanation for the observed 
brightness and polarization asymmetries, but (with only two clear- 
cut cases analysed in such detail so far) it is important to test it 
by looking at more objects. We r eviewe d the rest of the B2 low- 
luminosity source sample (Par ma et al. Ill987) to see if any other 
data su pport (or contest) the interpretation given here. Laing et al.l 
(I1999I) found that the source B2 0844+31 also has both a small jet 
to counter-jet width ratio and a high intensity ratio I -J Icy Unfor- 
tunately, there is no imaging for that source of the high quality we 
now have for 0206+35 and 0755+37 so we cannot test models of its 
asymmetries at the same level of detail. Nor can we classify its large 
scale structure definitively as Tobed' or 'plumed' : deeper imaging 
sensitive to its most extended structure is needed. Although lack of 
high-quality imaging precludes us from finding other good exam- 
ples of these phenomena in the B2 sample, we note that there are 
no clear counter-examples - either of sources in which the brighter 
jet appears to be wider than the counter-jet at low intensity lev- 
els, or of a large jet/counter-jet width asymmetry in a source that 
lacks 'lobed' structure or with only a small jet/counter-jet intensity 
asymmetry at its base^ 

As noted by iLaing et all (l2006al . see their fig. 15), the jets 
in the lobed FRI source 3C296 show /j//cj < 1 at their edges. 
The emission there is faint, but the effect is consistently present 
in the flaring and outer regions. The transverse variations of lin- 



ear polarization are also very different in the two jets ("Laing et aP 
[20p6a, Figs 18g and h): the counter-jet shows a prominent parallel- 
field edge, whereas the main jet does not. The model described by 
iLaing et all 12006a'). while giving a good overall fit to the bright- 
ness and polarization distributions of 3C 296, was not consistent 
with the observation of /j //cj < 1 and did not fully reproduce the 
flat profile of p with transverse apparent field in the approaching jet. 
We have examined possible backflow models for 3C 296 and find 
that they are qualitatively inconsistent with the polarization distri- 
bution, although they can easily fit the edge sidedness ratios. The 
combination of sidedness ratio and polarization is more reminiscent 
of the predictions of the outflow model analysed in Appendix lAl 

Emission from backflow such as that modelled here would be 
hard to recognise in lobed FRI sources whose jets are close to the 
plane of the sky. The backflow emission in such sources would be 
almost indistinguishable from faint outer edges of their jets and 
only unusually precise spectral index measurements could distin- 
guish it from low level brightness enhancements of the lobes near 
the jets. 

6.2 Should we expect backflows in FR I sources? 

Light jets propagating into dense media can be expected to termi- 
nate in one of two ways. They may decelerate and transition into 
'plumes' or 'tails' that are deflected away from the AGN by exter- 
nal pressure gradients or by winds in the IGM. Alternatively, they 
may deflect before reaching a contact discontinuity with the denser 
external medium, thus accumulating a 'cocoon' around the outflow. 
The first process is thought to underlie the formation of plumed or 
tailed FRI radio sources such as 3C 31 while the second is thought 
to form the 'classical double' lobed radio sources such as Cygnus A 
and is often associated with FR II morphology. Lobes in the (gener- 
ally more luminous) FRII sources also frequently contain discrete 
radio 'hot spots' that are identified with strong shocks where well- 
collimated (supersonic) outflows are slowed and begin to supply 
lobe material. There is no reason to suppose, however, that discrete 
hot spot formation is a requisite for cocoon (or radio lobe) forma- 
tion - momentum balance alone requires the deflection of the light 
outflow if it cannot escape along its initial path owing to devel- 
opment of a high pressure region downstream. Cocoons without 
hot spots are indeed seen in simulations of relativistic jets which 
are much li ghter than their surroundings (IPerucho & Martfl&007l : 
iRossi et aPboO S), in which the jets flows are transonic where they 
terminate. 

The majority of FRI sources form radio lobes whose detailed 
morphologies, spectral characteristics and polarization properties 
strongly resemble those of higher-power FR II lobes (iParma et al] 
ll996U1999l : lLaing et aLl20 1D. Their lobes have sharp outer bright- 
ness gradients, circumferential magnetic fields, and spectral indices 
that steepen towards the centre of the source on the largest angular 
scales - but without hot spots. Furthermore, outflows in lobed FRI 
source s can deflec t throu gh large angles without losing their iden- 
tities: iLaing et a l|) found regions where emission with jet- 
like spectral index ?^0.6 had displaced steeper- spectrum emission 
within FRI lobes. These results suggest that ongoing large-scale 
flow is present in these lobes well beyond the clearly recognisable 
jets. 

There is therefore both theoretical and observational support 
for supposing that jet outflows containing relativistic particles and 
magnetic fields may be redirected through large angles in lobed 
FRI sources. The additional ingredient suggested by our modelling 
of 0206+35 and 0755+37 is that a component of such an outflow 
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Figure 16. The fractional magnetic field components for the three sources, (a), (d) toroidal, {B^ / B'^)^/'^; (b), (e) longitudinal, {Bf / B'^y/'^; (c), (f) radial, 
(52/^2)1/2, (a) _ (c) 0206+35, (d) - (f) 0755+37. 



in an FR I source can return to the vicinity of the AGN as mildly 
relativistic backflow. As we noted in the introduction to this paper, 
this idea is supported by the presence of backflow with /3 > 0.2 
around the jets in some numerical simulations of the propagation 
of light, relativistic jets. The simulation by Perucho & Marti ( 2007|) 
used ini tial conditions for the jet derived from our FRI source 
models lLaing&Bridlell2002al lbl) and realistic density and pres- 
sure gradients in the s urrounding galactic and group atmosphere 
(iHardcastle et al.l l2002V In particular, the velocity at injection was 
P = 0.87 and the initial density contrast (the ratio of the density 
of the jet to that of its surroundings) was rj = 10~^. Although the 
jet had propagated only ^15 kpc by the end of the simulation, the 
structure already resembled a lobed FRI source of the type dis- 
cussed here, with a cocoon of backflowing, mixed jet and external 
plasma surrounding the jet. The jet was transonic at its termina- 
tion, so no hot spot was formed. Typical backflow velocities in the 
cocoon were /3 ^ 0.15, with values reaching [3 ^ 0.4 close to 
the nucleus. The use of an open boundary condition in the sym- 
metry plane at the base of the jet can cause the backflow s peed to 
be ove r-estimated dSaxton et al.ll2002h . although Peruc ho & Martil 
(l2007h argued that this effect was small in their simulation be- 
cause the flow through the open boundary was negligible. One other 
possible concern is that the simulation by Perucho & Marti (2007) 
was axisymmetric: the speed and extent of fast backflow appear to 
be smaller in some fully three-dimensio nal simulation s compared 
with the equivalent axisymmetric cases jNormanll 19961 : IaIov et aP 
11999'). We note, however, that the comparison may not be relevant 
to lobed FRI sources because the density contrast, 77 = 0.01, was 
much higher in these two examples, leading to cocoons which were 
far longer and thinner than those observed. The three-dimen siona l 
simulation of a relativistic jet with 77 = 10~^ bv lRossi et al.l ( l2008h 
indeed showed fast backflow with /3 ^ 0.4, despite the use of sym- 
metric boundary conditions at the jet inlet. The initial conditions 
(jet Lorentz factor F = 10) and the assumption of a uniform ex- 



ternal density are probably more appropriate to smaller physical 
scales than we consider here, however. Thus, alth ough the assump- 
tions a nd ii iitial conditions of the simulations by IPerucho & Martil 
(l2007h and lRossi et al ] 12008^ are not realistic enough to permit a 
quantitative comparison with our results, they do suggest that the 
idea of fast backflow is a reasonable one provided that the density 
contrast is very small ( < 10~^). 

The simulations discussed above are entirely hydrodynamic. 
We also note that backflow is an expected ingredient of models of 
magnetic hoop stress collimation of current-carrying jets because 
such models must provide a return current path - although it is 
unclear that such return paths need be as close to the jet outflow 
boundary as the backflow we have described here. 



7 SUMMARY AND FURTHER WORK 
7.1 Summary 

We have shown that many aspects of the intensity and linear po- 
larization distributions over the inner jets and counter-jets in the 
lobed FRI radio sources 0206+35 and 0755+37 are accounted for 
by an intrinsically symmetrical decelerating relativistic jet model 
that includes (mildly) relativistic backflow around both jets. 

We have estimated properties of this backflow subject to the 
simplifying assumptions that it is symmetrical across the AGN, ax- 
isymmetric, and that its streamlines are similar in shape to those 
of the outflow. Although these assumptions are likely to be too 
simple a priori we nevertheless find that the quality of the IQU 
fits obtained with the models including such symmetric backflow 
is similar to that o btained with pure dece l erating outflow models 
of other FRI jets ("Lains & Bridle 2002a!; ICanvin & La"in3 12004 
Canvin et al. 2005; Laing et al. 2006a). Furthermore, the outflow 
components of the models we have fitted to 0206+35 and 0755+37 
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are quite similar to those obtained for other FRI sources. The ad- 
dition of backflow to the models therefore suffices to explain the 
otherwise anomalous jet/counter-jet asymmetries of both sources 
and eliminates the need to invoke ad hoc environmental (or other 
intrinsic) side-to- side asymmetries. 

The salient features of backflow inferred from this procedure 
are as follows. 

(i) The backflow velocities are mildly relativistic, in the range 
0.05 </3< 0.35 (Fig.[Tl. 

(ii) The backflows are approximately symmetric around the out- 
flows and their radio emission comes from a hollow cone surround- 
ing the jet axis with additional half-opening angles ^ 8°. 

(iii) They can be traced to considerable distances from the AGN 
(at least 15 kpc for 0206-^35 and 50 kpc for 0755-^37) but the emis- 
sion close to the ends of the jets in both sources is chaotic, and it is 
not clear where the backflows begin. 

(iv) They do not emit synchrotron radiation all the way in to the 
AGN (Fig. US]). 

(v) The backflows emit with a spectral index a ^ 0.55 (Sec- 
tion [53]). This spectral index is lower than that of the nearby lobes 
and comparable with those of the outflows. 

(vi) Their magnetic fields are mostly toroidal and their emissiv- 
ities decrease with distance roug hly as r-^-^ (Figs [B] and [16]). 

These are the only two lobed FRI sources for which we 
have deep enough imaging and polarimetry to reveal the 'two- 
component' aspect of the jets and counter-jets that motivated this 
study. The generality of our results could thus be called into ques- 
tion by a single new example of an FRI source with strong jet- 
width asymmetries in which either (a) the axis is inferred to be 
close to the plane of the sky or (b) the apparently wider features are 
associated with the brighter jet. With only two examples of possi- 
ble backflow features we also cannot address whether all lobed FR I 
sources might contain backflow or (conversely) whether backflow 
exists only in lobed sources. 

The interpretation including backflow will continue to be 
preferable to any involving intrinsic side- to- side width differences 
if further studies find the apparently wider features only on the 
counter-jet side, and only in lobed sources for which inclination 
indicators suggest that the jets are at moderately large angles to the 
plane of the sky. 



7.2 Open questions and further work 

Our observations and models give no clue about the ultimate fate 
of the backflow or how it may interact with the outflow, but they 
raise a number of questions which could be addressed by deeper, 
higher-resolution observations of 0206-1-35 and 0755-1-37. 

(i) Where does the backflow originate? Does it start in a high- 
pressure region at the end of the outflow? 

(ii) Does the backflow shield the jet from entrainment or inter- 
action with the lobe plasma? 

(iii) Does the presence of the backflow perturb the jet structure 
in any way? 

(iv) Where does the backflow ultimately go: sideways or even 
closer to the AGN? 

(v) Why does the backflow radiate strongly where it does and 
stop radiating close to the AGN? 

(vi) Can the backflow really be faster than the asymptotic ve- 



locity of the outflow, as appears at first sight to be the case in 
0755-h3# 

Additional questions which could be answered by observa- 
tions of a sample of FR I sources include the following. 

(i) Do jets in other FRI sources with large jet- width asym- 
metries also have the two-component jet and counter-jet struc- 
ture found here in 0206-1-35 and 0755-1-37 (i.e. a strongly centrally 
brightened peaked main jet and centre-darkened counter-jet near 
the axis, and counter-jet emission consistently brighter than that of 
the main jet further from the axis)? 

(ii) Does the counter-jet/jet width asymmetry indeed correlate 
well with orientation indicators - counter-jet/jet intensity ratios 
and normalized core power - as expected in a relativistic backflow 
model of this asymmetry? If so, the tightness of the correlation with 
orientation indicators could be used to constrain the intrinsic sym- 
metry of the backflow. 

(iii) Does the width asymmetry indeed occur only in lobed 
FRI's? It will be important to obtain images which are sufficiently 
sensitive to extended structure to detect faint lobe emission in any 
sources whose structural classification is dubious. 

The high sensitivity and resolution of the imaging needed to ad- 
dress all of these issues and to test backflow models of the type we 
have proposed will require the use of the Jansky (Expanded) Very 
Large Array and e-MERLIN. 

Given the similarities between the extended emission in FRII 
and lobed FR I sources, it would also be interesting to search for ev- 
idence of backflow in the former class. The jets in FRII sources are 
usually much narrower than those we have imaged in the present 
study and are thought to be highly supersonic where they terminate 
in compact hot spots. Backflow is predicted by simulations of FR II 
dynamics, but it is unclear how its properties might depend on den- 
sity contrast, Mach number, magnetization and source age. It may 
be that observations of FR II sources without prominent hot spots 
will offer the best chance of detecting backflows. Counter-jets in 
FR II sources are faint and difficult to distinguish from filamentary 
lobe emission, so identification of any backflow component may be 
even more challenging than in FRI's. 

Three-dimensional simulations of very light, relativistic jets 
propagating in realistic external density and pressure distributions 
would be extremely valuable in understanding the backflow phe- 
nomenon in FRI sources. To be realistic, such simulations should 
be bipolar, with initial density contrasts ?^10~^. The effects of 
magnetic fields (ordered or disordered) on the flow also remain to 
be investigated. 



ACKNOWLEDGEMENTS 

The National Radio Astronomy Observatory is a facility of the Na- 
tional Science Foundation operated by Associated Universities, Inc. 
under co-operative agreement with the National Science Founda- 
tion. We are grateful to the referee for a very careful reading of 
the paper. RAL would also like to thank Alan and Mary Bridle for 
hospitality. 



^ The asymptotic outflow velocity is poorly constrained (Section [22), so 
this difference may not be real. 
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APPENDIX A: PURE OUTFLOW MODELS WITH 

/j//cj < 1 

It is possible under some circumstances for the ratio /j//cj (ap- 
proaching/receding) to be significantly less than unity close to the 



edges of the brightness distribution even for a symmetrical outflow. 
This might easily be mistaken for the effects of a backflowing com- 
ponent. We argue in this Appendix that the effect is quite likely to 
be observed in FRI jets, but that it is qualitatively inconsistent with 
the observations of 0206-1-35 and 0755-1-37 (particularl y in linear 
polarization). 3C 296 (modelled as a pure outflow by ILaing et al.l 
2006a) may show this effect at low brightness levels. 

/j//cj can become < 1 because of the effect of aberration 
on anisotropic rest-frame emission. If this acts in such a way that 
the magnetic field is nearly parallel to the line of sight in the rest 
frame, then the synchrotron emissivity can become very low. If this 
happens in the approaching jet but not the receding one, then the 
effect may be larger than that of Doppler boosting. A symmetrical 
pair of jets with purely toroidal fields can show this effect for some 
ranges of velocity. If the condition /3 = cos is satisfied at the edge 
of the approaching jet, then the field will be exactly parallel to the 
line of sight in the rest frame, so the synchrotron emissivity will be 
exactly zero. The condition can never be satisfied in the receding jet 
(except in the trivial case of a side-on source with zero velocity), 
so the sidedness ratio is also zero. Close to the jet edge or if the 
velocity condition is approximately satisfied, the sideness ratio can 
still be significantly less than unity. 

In order to demonstrate the effect, we consider a simple 
model with symmetrical, cylindrical, constant- velocity jets contain- 
ing purely toroidal fields. We also take the magnetic field and ra- 
diating particle density to be constant and assume a = 1 so that 
the calculated emi ssion p r ofiles are analytical, as given in the non- 
relativistic case bv lLaind(ll98ll) . Suppose that x is a coordinate in 
the plane of the sky perpendicular to the projected jet axis and nor- 
malized by the jet radius. Then the transverse profiles of sidedness 
and Q// are given by: 
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where the Doppler factors for the approaching and receding j ets are 

Dj = [r(l -/3cos6')]"^ (A4) 
Dei = [r(l + /3cos6')]-' (A5) 



We show some example profiles in Fig. lAll We have established 
that the magnetic-field structures of FRI jets tend to be toroidally 
dominated at large distances from the nucleus and Fig. lAll shows 
that /j//cj < 1 at the edges for plausible velocities and angles 
to the line of sight, so it would not be surprising to see this ef- 
fect in some sources. An inevitable corollary, however, is that the 
degree of polarization at the edges of the jets must be high, with 
the apparent magnetic field transverse to the jet axis. The rea- 
son is that the toroidal field loops are seen close to edge-on in 
the rest frame: in particular, we should not observe the transition 
from transverse apparent field on-axis to longitudinal at the edges. 
The main jet transverse profiles of Q/I for 0206-1-35 and 0755-1-37 
(Figs [8] and O indicate that the apparent field is primarily lon- 
gitudinal (Q/I < 0) at the edges and certainly inconsistent with 
the predicted Q/I ^ +0.7. The sources must also have ^<40° 
in order to produce the large values of /j//cj observed for their 
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Figure Al. Transverse profiles of sidedness ratio and Qjl for cylindri- 
cal model jets containing purely toroidal fields (equations IAll - lA5l . (a) 
and (d): sidedness ratio I-Jicy (b) and (e): (main jet), (c) and (f): 

Qcj/^cj (counter-jet), (a) - (c): Fixed angle to the line of sight, 6 = 60°. 
The velocities are /5 = 0.1 (dots), 0.25 (short dash), 0.5 (full) and 0.75 
(long dash), (d) - (f): Fixed velocity /3 = 0.1 The angles to the line of sight 
are = 40° (dots), 60° (short dash), 80° (short dash) and 84.23° (full). 
In all of the panels, the full lines represent the case (3 = cos 6, for which 
toroidal field loops are seen edge-on in the rest frame in the main jet, so 
Qj//j =+po = +3/4. 



jet bases. This in turn requires high edge velocities to satisfy the 
condition (3 ^ cosO, giving a very narrow edge with /j//cj < 1. 
Detailed modelling confirms that the predicted brightness and po- 
larization distributions are quite unlike those observed in 0206+35 
and 0755+37. 

The sidedness and Q/I profi les are, however, qualitatively 
similar to those observed in 3C 296 (iLaing et al.l2006ah . except that 
the observed value of Q/I for the main jet of 3C 296 is ^0.3, com- 
pared with the predicted 0.7 for a pure toroidal field. Detailed mod- 
elling confirms that a simple field configuration of this type cannot 
simultaneously fit the sidedness ratio and polarization, but the edge 
emission is very faint so contamination by lobe emission may be 
significant: deeper observations are needed in order to separate jet 
and lobe emission unambiguously. 



